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Abstract 

We present a mathematical study of the order book as a multidimensional continuous- 
time Markov chain where the order flow is modeled by independent Poisson processes. Our 
aim is to bridge the gap between the microscopic description of price formation (agent- 
based modeling), and the Stochastic Differential Equations approach used classically to 
describe price evolution at macroscopic time scales. To do this, we rely on the theory 
of infinitesimal generators and Foster-Lyapunov stability criteria for Markov chains. We 
motivate our approach using an elementary example where the spread is kept constant 
( "perfect market making" ) . Then we compute the infinitesimal generator associated with 
the order book in a general setting, and link the price dynamics to the instantaneous 
state of the order book. In the last section, we prove that the order book is ergodic— 
in particular it has a stationary distribution — that it converges to its stationary state 
exponentially fast, and that the large-scale limit of the price process is a Brownian motion. 

Keywords: Limit order book; agent-based modeling; order flow; bid-ask spread; Markov 
chain; stochastic stability; FCLT; geometric mixing. 

1 Introduction and Background 

The emergence of electronic trading as a major means of trading financial assets makes the 
study of the order book central to understanding the mechanisms of price formation. In order- 
driven markets, buy and sell orders are matched continuously subject to price and time priority. 
The order book is the list of all buy and sell limit orders, with their corresponding price and 
size, at a given instant of time. Essentially, three types of orders can be submitted: 

• Limit order: Specify a price (also called "quote") at which one is willing to buy or sell a 
certain number of shares; 

• Market order: Immediately buy or sell a certain number of shares at the best available 
opposite quote; 
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• Cancellation order: Cancel an existing limit order. 



In the econophysics literature, "agents" who submit exclusively limit orders are referred to as 
liquidity providers. Those who submit market orders are referred to as liquidity takers. 

Limit orders are stored in the order book until they are either executed against an incoming 
market order or canceled. The ask price P A (or simply the ask) is the price of the best (i.e. 
lowest) limit sell order. The bid price P B is the price of the best (i.e. highest) limit buy order. 
The gap between the bid and the ask 

S: =p A -P B } (1) 

is always positive and is called the spread. Prices are not continuous, but rather have a discrete 
resolution AP, the tick, which represents the smallest quantity by which they can change. We 
define the mid-price as the average between the bid and the ask 

pA + pB 

* = — (2) 

The price dynamics is the result of the interplay between the incoming order flow and the 
order book [2J. Figure [I] is a schematic illustration of this process [I]. Note that we chose to 
represent quantities on the bid side of the book by non-positive numbers. 

Although in reality orders can have any size, we shall assume throughout this paper that 
all orders have a fixed unit size q. This assumption is convenient to carry out our analysis and 
is, for now, of secondary importance to the problem we are interested in. 



2 An Elementary Approximation: Perfect Market Mak- 
ing 

We start with the simplest agent-based market model: 

• The order book starts in a full state: All limits above P A (0) and below P B (0) are filled 
with one limit order of unit size q. The spread starts equal to 1 tick; 

• The flow of market orders is modeled by two independent Poisson processes M + (t) (buy 
orders) and M~(t) (sell orders) with constant arrival rates (or intensities) A + and A~; 

• There is one liquidity provider, who reacts immediately after a market order arrives so 
as to maintain the spread constantly equal to 1 tick. He places a limit order on the same 
side as the market order (i.e. a buy limit order after a buy market order and vice versa) 
with probability u and on the opposite side with probability 1 — u. 

The mid-price dynamics can be written in the following form 

dP{t) = AP(dM + (t) - dM~(t))Z, (3) 
where Z is a Bernoulli random variable 

Z = with probability (1 — u), 
Z = 1 with probability u. 
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(1 ) initial state (2) liquidity is taken (3) wide spread 




(4) liquidity returns (5) liquidity returns (6) final state 
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Figure 1: Order book schematic illustration: a buy market order arrives and removes liquidity 
from the ask side, then sell limit orders are submitted and liquidity is restored. 
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The infinitesimal generatoiQ associated with this dynamics is 

Cf(P) = u (X + (f(P + AP) -f) + \-(f(P - AP) - /)) . (6) 

It is well known that a continuous limit is obtained under suitable assumptions on the intensity 
and tick size. Noting that ^ can be rewritten as 

Cf(P) = 1 -u(X+ + X-)(AP? f{P + AP) - 2f + fiP - AP) 



2 v yv ' (AP) 
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+ m (a + -a W (P + AP >-/ (P " AP) ' <*> 



and under the following assumptions 

u(A+ + \-){AP) 2 — >a 2 as AP ^ 0, 

u(X + - \-)AP^/2 asAP^O, ((N) 

the generator converges to the classical diffusion operator 

a 2 d 2 f df_ 

corresponding to a Brownian motion with drift. This simple case is worked out as an example 
of the type of limit theorems that we will be interested in in the sequel. One should also note 
that a more classical approach using the Functional Central limit Theorem (FCLT) as in pQ or 
[T2] yields similar results ; For given fixed values of A + , A~ and AP, the rescaled-centered price 
process 

P(nt) — nfit 



ma 

converges as n — > oo, to a standard Brownian motion (B(t)) where 



(10) 



a = APy/(X+ + X-)u, 
fi = AP(A + - \~)u. 

Let us mention that one can easily achieve more complex diffusive limits such as a local volatility 
model by imposing that the limit is a function of P and t 

u(X + + \-)(AP) 2 ^a 2 (pt), 

u(X + - X-)AP -)• fi(P,t). 1 ' 

This would be the case if the original intensities are functions of P and t themselves. 



1 The infinitesimal generator of a time-homogeneous Markov process (X.(t))t>o is the operator £, if exists, 
defined to act on sufficiently regular functions / : K n — >• K, by 

, ,. E[/(X(t))|X(0)=x]-/(x) 
£/(x) := hm 1 . (5) 

It provides an analytical tool to study (X(t)) [6]. 
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3 Order Book Dynamics 



Model setup: Poissonian arrivals, reference frame and boundary con- 
ditions 

We now consider the dynamics of a general order book under a Poisson type assumption for 
the arrival of new market orders, limit orders and cancellations. We shall assume that each 
side of the order book is fully described by a finite number of limits N, ranging from 1 to iV 
ticks away from the best available opposite quote. We will use the notation 

X(t) := (a(t); b(*)) := (ai(t), • • • , a N (t); b±{t), . . . , b N {t)) , (13) 

where a:=(ai,...,ajv) designates the ask side of the order book and a% the number of shares 
available at price level i (i.e. i ticks away from the best opposite quote), and b := (pi, . . . ,6jv) 
designates the bid side of the book. By doing so, we adopt the representation described e.g. 
in [3] or but depart slightly from it by adopting a finite moving frame, as we think it is 

realistic and more convenient when scaling in tick size will be addressed. 
Let us now recall the events that may happen: 

• arrival of a new market order; 

• arrival of a new limit order; 

• cancellation of an already existing limit order. 

These events are described by independent Poisson processes: 

• M ± {t): arrival of new market order, with intensity A M+ I(a ^ 0) and A M I(b ^ 0); 

• Lf(t): arrival of a limit order at level i, with intensity Xf ; 

• Cf(t)\ cancellation of a limit order at level i, with intensity Xf + — and Xf — — . 

q is the size of any new incoming order, and the superscript "+" (respectively "— ") refers to 
the ask (respectively bid) side of the book. Note that the intensity of the cancellation process 
at level i is proportional to the available quantity at that level. That is to say, each order at 
level % has a lifetime drawn from an exponential distribution with intensity Xf . Note also that 
buy limit orders L^(t) arrive below the ask price P A (t), and sell limit orders Lf(t) arrive above 
the bid price P B {t). 

We impose constant boundary conditions outside the moving frame of size 2N: Every time 
the moving frame leaves a price level, the number of shares at that level is set to fl m (or b^ 
depending on the side of the book). Our choice of a finite moving frame and constant]^] boundary 
conditions has three motivations. Firstly, it assures that the order book does not empty and 
that P A , P B are always well defined. Secondly, it keeps the spread S and the increments of 



2 See also [5] for an interesting discussion. 

3 Actually, taking for and \boo \ independent positive random variables would not change much our analysis. 
We take constants for simplicity. 
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Figure 2: Order book dynamics: in this example, N = 9, q = 1, a c 
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-4. The shape of 



the order book is such that a(t) = (0, 0, 0, 0, 1, 3, 5, 4, 2) and b(t) = (0, 0, 0, 0, -1, 0, -4, -5, -3). 
The spread S(t) = 5 ticks. Assume that at time t' > t a sell market order dM~(t') arrives, 
then a(t') = (0,0,0,0,0,0,1,3,5), b(t') = (0,0,0,0,0,0,-4,-5,-3) and S(t') = 7. Assume 
instead that at t' > t a buy limit order dL^(t') arrives one tick away from the best opposite 
quote, then a(t') = (1,3,5,4,2,4,4,4,4), b(t') = (-1,0,0,0,-1,0,-4,-5,-3) and S(t') = 1. 



P A , P B and P = p ^ bounded — This will be important when addressing the diffusive limit 
of the price. Thirdly, it makes the model Markovian as we do not keep track of the price 
levels that have been visited (then left) by the moving frame at some prior time. Figure [2] is a 
representation of the order book using the above notations. 
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Evolution of the order book 

We can now write the following coupled SDEs for the quantities of outstanding limit orders in 
each side of the order bookEl 

daft) = -[q-J2 ak ) dM+ ( t ) + i dL t(t) ~ QdCf{t) 
V fc=i / + 

N N 

+ (J M "(a) - a)idM-(t) + (J L <(*) ~ a)i^-(t) + £ (J cr (a) - a)^Cr(t), 

i=l i=l 

tfti(0 = ( 9 - X) |6fc| ) dM~(*) - gdLr(t) + gdCr(t) 

V fe=i / + 

AT TV 

+ (J M+ (b) - b)idM+(t) + X (^( b ) - *>)idLt(t) + X ( jCl+ ( b ) " b)^C+(t), 

i=l i=l 

(14) 

where the J's are s/w/i operators corresponding to the renumbering of the ask side following 
an event affecting the bid side of the book and vice versa. For instance the shift operator 
corresponding to the arrival of a sell market order dM~(t) of size q iaj 

J M ~ (a) = I 0, 0, . . . , , ax, a 2 , . . . , aN-k J , (15) 

\ k times / 

with 

v 

k := ini{p : Ifyl > <?} ~ inf {P : \K\ > °l- 

i=i 

Similar expressions can be derived for the other events affecting the order book. 

In the next sections, we will study some general properties of such models, starting with 
the generator associated with this 2iV-dimensional continuous-time Markov chain. 



4 Rcmember that, by convention, the bi's are non-positive. 

5 For notational simplicity, we write J M (a) instead of J M (a; b) etc. for the shift operators. 
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4 Infinitesimal Generator 



Let us now work out the infinitesimal generator associated with the jump process described 
above. We have 

Cf (a; b) = X M+ (f {[a, -(q- A(i - l)) + ] + ; J M+ (b)) - /) 

N 

+ E A * L+ (/ (ai + q;J L Hb))-f) 

i=i 

+ f>? + % (a,-,;J^ + (b))-/) 



<7 

i=i 

\M~ I f I jM 



+ A (/ [6, + (g - S(i - l)) + ]_ ) - / 

AT 

+ E A H/ (a) ;6i -?)-/) 



i=i 

N 



1=1 ^ 



(17) 



where, to ease the notations, we note /(o^; b) instead of f(a±, . . . , a^, . . . , a^; b) etc. and 

x + := max(x, 0), x_ := min(x,0),x G K. (18) 

The operator above, although cumbersome to put in writing, is simple to decipher: a series of 
standard difference operators corresponding to the "deposition-evaporation" of orders at each 
limit, combined with the shift operators expressing the moves in the best limits and therefore, 
in the origins of the frames for the two sides of the order book. Note the coupling of the two 
sides: the shifts on the a's depend on the 6's, and vice versa. More precisely the shifts depend 
on the profile of the order book on the other side, namely the cumulative depth up to level % 
defined by 



k=l 



(19) 



k=i 



and the generalized inverse functions thereof 

p 

A-\q') :=mf{p:J2^>q'}, 



3=1 
P 



(20) 



B-\q') :=inf{p:^|6,| > q'}, 



where q' designates a certain quantity of shares. Note that a more rigorous notation would be 

A(i,SL(t)) and A _1 (g',a(t)) 
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for the depth and inverse depth functions respectively. We drop the dependence on the last 
variable as it is clear from the context. 

Remark 4.1 The index corresponding to the best opposite quote equals the spread S in ticks, 
that is 



% A := A-\0) = M{p : E a i > °) 



v 



AP 



i B ■= B-\0) = M{p : N > 0} = 



S 
AP 



:= t s = i A - 



(21) 



5 Price Dynamics 

We now focus on the dynamics of the best ask and bid prices, denoted by P A {t) and P B {t). 
One can easily see that they satisfy the following SDEs 

' dP A (t) = APKA-^q) - A-\0))dM + (t) 

N 

- E " *)+ dL t(t) + (A~\ q ) - A-\0))dC+(t)}, 



i=l 



dP B (t) = -AP^" 1 ^) - B-\<S))dM-(t) 



(22) 



A' 



which describe the various events that affect them: change due to a market order, change due 
to limit orders inside the spread, and change due to the cancellation of a limit order at the 
best price. One can summarize these two equations in order to highlight, in a more traditional 
fashion, the respective dynamics of the mid-price and the spread 



dP(t) 



AP 

N 



[{A-\q) - A-\<d))dM + (t) - (B-\q) - B' 1 (0))dM- (t) 



N 



- E - + E - 

+ (A-\q) - A-\0))dCt A (t) - (B-\q) - B-\0))dCr(t)] , 



(23) 



dS(t) = AP [(A'^q) - A-\<d))dM + {t) + (B-\q) - B-\G))dM-(t) 

N N 

- E (^(O) - i) + dL+(t) - E (B~\0) - i) + dL-(t) 

i=i i=i 
+ (A-\q) - A-\0))dC+(t) + (B-\q) - B-\0))dCr(t)] . 



(24) 



The equations above are interesting in that they relate in an explicit way the profile of the 
order book to the size of an increment of the mid-price or the spread, therefore linking the 
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price dynamics to the order flow. For instance the "infinitesimal" drifts of the mid-price and 
the spread, conditional on the shape of the order book at time t, are given by 



E[dP(t)|(a;b)] 



AP 

N 



(A-\q) - A-\Q))\ M - - {B-\q) - 5- i (0))A 



M~ 



N 



i=l 



i=l 



and 



E [dS(t)\( a ;b)} = AP I (A-\q) - A-\0))\ M+ + {B-\q) - B~\Q))\ 



N 



N 



- E (^P) - ^ x " + - E - ^ 



dt. 



i=i 

-\. i ,.\ i-[,n,nO+°»A i ft3-lf„\ n>-l/n\\ \C~ l^sl 



+ (^(g) - A-(0))A^^ + (B-\q) - 5- 1 (0))A^ f 



(25) 



cit. (26) 



6 Ergodicty and Diffusive Limit 

In this section, our interest lies in the following questions: 

1. Is the order book model defined above stable? 

2. What is the stochastic-process limit of the price at large time scales? 

The notions of "stability" and "large scale limit" will be made precise below. We first need 
some useful definitions from the theory of Markov chains and stochastic stability. Let (Q*)t>o 
be the Markov transition probability function of the order book at time t, that is 

Q*(x, E) := P [X(t) G E\X(0) = x] , t G R+, x G S, E C S, (27) 

where S C Z 2JV is the state space of the order book. We recall that a (aperiodic, irreducible) 
Markov process is ergodic if an invariant probability measure ir exists and 



lim ||Q'(x, 



7T 



0,Vx G 5, 



t— s>oo 



where ||.|| designates for a signed measure v the total variation normH defined as 



\u\\ := sup |^(/)| = sup v{E) — inf //(£"). 
/:|/|<i £eB(S) £eB(s) 



3 The convergence in total variation norm implies the more familiar pointwise convergence 

lim |Q t (x,y)- 7 r(y)| =0,x,y G 5. 

t—¥oo 



(2f 



(30) 



(29) 



Note that since the state space 5 is countable, one can formulate the results without the need of a "measure- 
theoretic" framework. We prefer to use this setting as it is more flexible, and can accommodate possible 
generalizations of our results. 
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In (30), B(S) is the Borel cr-field generated by S, and for a measurable function / on S, 
u(f):=J s fdu. 

V -uniform ergodicity. A Markov process is said V -uniformly ergodic if there exists a coer- 
civ^] function V > 1, an invariant distribution tt, and constants < r < 1, and R < oo such 
that 

||Q*(x,.)-7r(.)|| <&V(x),xG5,i ER+. (31) 

V— uniform ergodicity can be characterized in terms of the infinitesimal generator of the Markov 
process. Indeed, it is shown in [7J [8] that it is equivalent to the existence of a coercive function 
V (the "Lyapunov test function" ) such that 

CV(x) < -pV(x) + 7, (Geometric Drift condition.) (32) 



for some positive constants (3 and 7. (Theorems 6.1 and 7.1 in [8].) Intuitively, condition (J32| 
says that the larger V(X(£)) the stronger X is pulled back towards the center of the state space 
S. A similar drift condition is available for discrete-time Markov processes (X n )„ g N and reads 

W(x) < -/3K(x)+ 7 I c (x), (33) 

where T> is the drift operator 

W(x) := E[V(X n+1 ) - K(X n )|X n = x]. (34) 
and C C S a finite set. (Theorem 16.0.1 in [7j.) We refer to [7] for further details. 

Ergodicity of the order book and rate of convergence to the stationary 
state 

Of utmost interest is the behavior of the order book in its stationary state. We have the 
following result: 

Theorem 6.1 If Xg = mini<j<Ar{Af ± } > 0, then (X(t)) t > = (a(t); h(t)) t >o is an ergodic 
Markov process. In particular (X(i)) has a stationary distribution tt. Moreover, the rate of 
convergence of the order book to its stationary state is exponential. That is, there exist r < 1 
and R < 00 such that 

||Q*(x,.) -tt(.)|| <^(x),tGM + ,xG5. (35) 

Proof. Let 

Af Af 

V(x) := y (a; b) := ^ a 1 + ^ |^| + q (36) 



i=i i=i 



7 That is, a function such that V(x) — > 00 as |x| — > 00. 
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be the total number of shares in the book (+q shares). Using the expression of the infinitesimal 
generator (17) we have 



N 



N 



CV{x) < —(X M+ + \ M ~)q + J>f + Af) ? - $>fV + Af |6*|) 



N 



N 



L+, 



1*5 



L *'ig-i) + \U"\ 



i=l 



< -(A A/+ + \ M )q + (A^~ + A^)g - A^(x) 



+ N(A L ~a°° + A L+ |6°°|) 



where 



N 



A L± := V Af and A c := min {Af } > 0. 

' Ki<JV 



(37) 



(3* 



(39) 



i=l 



The first three terms in the right hand side of inequality (37) correspond respectively to the 
arrival of a market, limit or cancellation order — ignoring the effect of the shift operators. The 
last two terms are due to shifts occurring after the arrival of a limit order inside the spread. 
The terms due to shifts occurring after market or cancellation orders (which we do not put in 
the r.h.s. of (37)) are negative, hence the inequality. To obtain inequality (38), we used the 
fact that the spread is is bounded by + 1 — a consequence of the boundary conditions we 
impose — and hence (is — i)+ is bounded by N. 



The drift condition ( 38 ) can be rewritten as 

£V(x) < -/3V(x)+ 7 , 



(40) 



for some positive constants /3, 7. Inequality (40) together with theorem 7.1 in [H] let us assert 
that (X(i)) is ^-uniformly ergodic, hence (35). □ 

Corollary 6.1 The spread S(t) = v4 _1 (0, a(i))AP = 5(X(t)) has a well-defined stationary 
distribution — This is expected as by construction the spread is bounded by N + 1. 



The embedded Markov chain 

Let (X n ) denote the embedded Markov chain associated with (X(£)). In event time, the prob- 
abilities of each event are "normalized" by the quantity 

N N 

A(x) := A M+ + X M ~ + A L+ +A L ~+J2 Af * + E A f N- ( 41 ) 

i=l i=l 

For instance, the probability of a buy market order when the order book is in state x, is 

X M+ 

P["Buy market order at time n" |X n _i = x] := p M+ (x) = — -. (42) 
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The choice of the test function V(x) = J2i a i + J2i k + Q does not yield a geometric drift 
condition, and more care should be taken to obtain a suitable test function. Let z > 1 be a 
fixed real number and consider the function^ 

V(x) : = zE^+Eil 6 *! := ^ (x) . (43) 

We have 

Theorem 6.2 (X n ) is V -uniformly ergodic. Hence, there exist r2 < 1 and R% < oo such that 

||f/ n (x,.) <% 2 n y(x),neN,,x6 5. (44) 

where (U n ) ne ^ is the transition probability function o/(X„) ne N and v its stationary distribution. 
Proof. 

W(x) < ^— -Wx)) 
~~ A(x) v v 77 

^ A(x) v 



X M 



+ 2[ ) ai+jva °° +y:w i&ii+<? v(x)) 



A(x 

If we factor out V(x) = ^D^+Z)^ i n the r.h.s of (45), we get 

W(x) £ A"++A"- 



A C ~I6 I 

+ J ' J ' (^ a ' +E ' |fe ' hg -V(x)). (45) 



K(x) " A(x) 
A L ~ + A L 



A(x) ^ 11 

+ ^V^N ( ^- 1)t (46) 

where 

doo := max{aoo, |&oo|}. (47) 



3 To save notations, we always use the letter V for the test function. 
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Then 

W(x) < X^+X^ _ 

V{*) ~ \M+ + \M- +A L+ +A L- +X C^) K 

i A + A (~q+Ndoo _ 1 \ 

A M+ + A M ~ + A i+ + A L ~ + A£<£>(x) ; 

+ ^ = (.--1), (48) 

A M+ + A M ~ + A L+ + A L ~ + \ c (p(x) 

with the usual notations 

A^ := min Xf ± and A^ := max Af* . (49) 



Denote the r.h.s of (48) -B(x). Clearly 



lim B(x) = =U= J < 0, (50) 

hence there exists A > such that for xG<S and y(x) > A 

W(x) A c (2~ 9 -1) „ , N 

—V < — — — " := -P < 0. 51 

V(x) 2X C K ' 

Let C denote the finite set 

C = {xG5:^(x) = ^a 1 + ^5 i <4 (52) 

i i 

We have 

W(x) < -/3F(x)+ 7 I c (x), (53) 

with 

7 := maxW(x). (54) 
Therefore (X n ) n > is ^-uniformly ergodic, by theorem 16.0.1 in [7]. □ 

The case of non-proportional cancellation rates 

The proof above can be applied to the case where the cancellation rates are independent of 
the state of order book X'(i) — We shall denote the order book X'(£) in order to highlight that 
the assumption of proportional cancellation rates is relaxed. The probability of a cancellation 
dCf(t) in [t,t + 6t] is now 

P[Cf {t + St) - Cf{t) = l|X'(t) = x'] = xfSt + o(St), (55) 

where lim 5 ^ o{St)/St = 0. Since A = X M+ + X M ~ + A L+ + A L ~ + J2?=i A? + + Eii A ?~ does 
not depend on x', the analysis of the stability of the continuous-time process (X'(£)) and its 
discrete-time counterpart (XJJ are essentially the same. 
We have the following result: 
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Theorem 6.3 Set 



Under the condition 



N 



N 



A c± :=^Af andA L± :=$>f. 



X M+ + X M ~ + A°' T + A u ~ > (A L ^ +A L ')(1 + N 



c+ 



i=i 



.c 



L+ 



(X' n ) is V -uniformly ergodic. There exist r 3 < 1 and R3 < 00 such that 

||t/'"(x,.) -!/'(.) 1 1 < ^3,^ 6N,xG5. 
T/ie same is true for (X!(t)). 



(56) 



(57) 



(5f 



Proof. Let us prove the result for (XJJ. Inequality (46) is still valid by the same arguments, 
but this time the arrival rates are independent of x' 



W(x') 



\M+ 1 \Af- 



+ 
+ 



Set 

A Taylor expansion in e gives 

W(x) 



A 

A L+ + A L 

A v ~ 
A c+ +A C ~ 
A {2 

z =: 1 + e > 1. 



y q+Nd a 



1) 



1)- 



(59) 
(60) 



A- 



V(x) 



< 


(A M+ 


+ A M ~ 


+ 


(A i+ 


+ A L ~) 


+ 


(A c+ 


+ A C ~ 



For e > small enough, the sign of (|61|) is determined by the quantity 



(A L+ + A L ~)(l + N-) - (A°' T + A c '~). 



Ac 



Hence, if (57) holds 



W(x) < -0V(x) for some /3 > 0, 
and a geometric drift condition is obtained for X'. 



(61) 
(62) 

(63) 

□ 



If for concreteness we set q = 1 share, and all the arrival rates are constant and symmetric, 

(64) 



then condition (57) can be rewritten as 

X M + NX C > NX L (1 + Nd 



where N is the size of the order book and is the depth far away from the mid-price. Note 
that the above is a sufficient condition for (^-uniform) stability. 
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Large scale limit of the price process 

We are now able to answer the main question of this paper. Let us define the process e(t) G 
{1, . . . , 2(2N + 1)} which indicates the last event 

{M ± ,Lf ,Cf }; G {i,...,jv}, 

that has occurred before time t. 

Lemma 6.1 If we append e(t) to the order book (X(t)) ; we get a Markov process 

Y(f):=(X(f),e(t)) (65) 



which still satisfies the drift condition (32). 

Proof. Since e(t) takes its values in a finite set, the arguments of the previous sections are 
valid with minor modifications, and with the test functions 

V(y) := q + a 1 + \b l \ + e, (continuous-time setting) (66) 
V(y) := e Eai+EN+e . (discrete-time setting) (67) 
The ^/-uniform ergodicity of (Y(t)) and (Y n ) follows. □ 

Given the state X n _! of the order book at time n — 1 and the event e n , the price increment at 
time n can be determined. (See equation (23).) We define the sequence of random variables 

Vn := *(X n _ 1; e n ) := $(Y n , Y n _x), (68) 

as the price increment at time n. ^ is a deterministic function giving the elementary "price- 
impact" of event e n on the order book at state X n _i. Let \x be the stationary distribution of 
(Y n ), and M its transition probability function. We are interested in the random sums 

n n 
k=l k=l 

where 

rj k := Vk - EM = $ fc = $ fe - E M [$ fc ], (70) 
and the asymptotic behavior of the rescaled-centered price process 

pWf t ) : = 3^L, (71) 

as n goes to infinity. 
Theorem 6.4 The series 

oo 

a 2 = E fl [f ] + 2j2^[Wn] (72) 

n=i 

converges absolutely, and the rescaled-centered price process is a Brownian motion in the limit 
of n going to infinity. That is 

p(n)( t ) aB (t), (73) 

where (B(t)) is a standard Brownian motion. 
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Proof. The idea is to apply the Functional Central Limit Theorem for (stationary and ergodic) 
sequences of weakly dependent random variables with finite variance. First, we note that the 
variance of the price increments r] n is finite since it is bounded by JV+1. Second, the V^-uniform 
ergodicity of (Y n ) is equivalent to 



|M n (x,.) < Rp n v(*),n G N, 



(74) 



for some R < oo and p < 1. This implies thanks to theorem 16.1.5 in [7J that for any g 2 , h 2 < V, 
k, n G N, and any initial condition y 

\E y [g(Y k )h(Y n+k )]-E y [g(Y k )]E y [g(Y k )]\ < Rp n [l + p k V(y)}, (75) 

where E y [.] means E[.|Yo = y]. This in turn implies 

\E y [h(Y k )g{Y k+n )}\ < R lP n [l + p k V(y)} (76) 

for some R\ < oo, where h = h — E M [/i],y = g — E^g]. By taking the expectation over p on 
both sides of (76) and noting that E A1 [V A (Y )] is finite by theorem 14.3.7 in [7], we get 

|E^(Y fc )#(Y fc+n )]| < R 2 p n =: p(n),k,n G N. (77) 
Hence the stationary version of (Y n ) satisfies a geometric mixing condition, and in particular 

J>(n)<oo. (78) 



Theorems 19.2 and 19.3 in pQ on functions of mixing processeqj let us conclude that 

oo 

a 2 :=E /i [^] + 2^E M [r/ fJ n ] (79) 



n=l 



is well-defined — the series in (79) converges absolutely — and coincides with the asymptotic 
variance 



lim — E. 

n— >oo n 



k=l 



a' 



Moreover 



iO) 



where (B(t)) is a standard Brownian motion. The convergence in (81 ) happens in -D[0, oo), the 
space of M-valued cadlag functions, equipped with the Skorohod topology. □ 



Remark 6.1 Obviously, theorem \6.J\ is also true with non-proportional cancellation rates under 
condition (57). In this case the result holds both in event time and physical time. Indeed, let 
(N(t))teSL + denote a Poisson process with intensity A = A + A + ^f =1 A? . The price 



process in physical time (P c (t))tm+ can be linked to the price in event time (P n )n&N by 



Pc(t) = P. 



N(t)- 



(82) 



9 See also theorem 4.4.1 in [12] and discussion therein. 
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Then 



aB(t) as in theorem 6.4, (83) 



1 ■ N(u) «->qo 1 

ana since — ^ — > 1 a.s., 



Remark 6.2 In i/ie /arge sca/e Zimit, £/ie mid-price P, the ask price P A — P + f , and i/ie 6id 
price P B = P- I converge to the same process (o~B(t)). 

Numerical illustration 



Figures [5}|8] are obtained by numerical simulation of the order book. We note in particular 
the asymptotic normality of price increments, the fast decay of autocorrelation and the linear 
scaling of variance with time, in accordance with the theoretical analysis. 

7 Conclusion and Prospects 

This paper provides a simple Markovian framework for order book modeling, in which elemen- 
tary changes in the price and spread processes are explicitly linked to the instantaneous shape 
of the order book and the order flow parameters. Two basic properties were investigated: the 
ergodicity of the order book and the large scale limit of the price process. The first property, 
which we answered positively, is desirable in that it assures the stability of the order book in the 
long run. The scaling limit of the price process is, as anticipated, a Brownian motion. A key 
ingredient in this result is the convergence of the order book to its stationary state at an expo- 
nential rate, a property equivalent to a geometric mixing condition satisfied by the stationary 
version of the order book. This short memory effect, plus a constraint on the variance of price 
increments guarantee a diffusive limit at large time scales. We hope that our approach offers a 
plausible microscopic justification to the much celebrated Bachelier model of asset prices. 

We conclude with a final remark regarding two possible extensions: The assumption of a 
finite order book size — and hence a bounded spread — may seem artificial, and one can seek 
more general stability conditions for an order book model in which the spread is unbounded 
a priori. In addition, richer price dynamics (heavy tailed return distributions, long memory, 
more realistic spread distribution etc.) can be achieved with more complex assumptions on the 
order flow (e.g. feedback loops [ID] , or mutually exciting arrival rates [9]). These extensions 
may, however, render the model less amenable to mathematical analysis, and we leave the 
investigation of such interesting (but difficult) questions for future research. 
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Figure 3: Average depth profile. 
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Figure 4: Histogram of the spread. 
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Figure 7: Autocorrelation of price increments. 
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Figure 8: Variance in event time. The dashed line is a linear fit. 
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